Mon. Not. R. Astron. Soc. 000, [T]-?? (2008) Printed 8 November 2012 (MN WF^K. style file v2.2) 



The "universal" radio/X-ray flux correlation : the case 
study of the black hole GX 339-4. 

g: S. Corbel 1 ' 2 * , M. Coriat 3 ' 1 , C. Brocksopp 4 , A.K. Tzioumis 5 , R. P. Fender 3 , 
c5: J.A. Tomsick 6 , M. M. Buxton 7 , C. D. Bailyn 7 

1 Laboratoire AIM (CEA/IRFU - CNRS/INSU - Universite Paris Diderot), CEA DSM/IRFU/SAp, F-91191 Gif-sur-Yvette, France 
2 Institut Universitaire de France, 75005 Paris, France. 

3 School of Physics and Astronomy, University of Southampton, Highfield, Southampton, SO 17 IB J, UK. 
4 Mullard Space Science Laboratory, University College London, Holmbury St Mary, Dorking, Surrey RH5 6NT, UK. 
5 Australia Telescope National Facility, CSIRO, P.O. Box 76, Epping NSW 1710, Australia. 
6 Space Sciences Laboratory, 7 Gauss Way, University of California, Berkeley, CA 94720-7450, USA. 
7 Astronomy Department, Yale University, P.O. Box 208101, New Haven, CT 06520-8101, USA. 



> 

o 



> 
o 



- 1—^ 

X 



November 2012 



ABSTRACT 



w 
X 

6 

H 

1 The existing radio and X-ray flux correlation for Galactic black holes in the hard 

, and quiescent states relies on a sample which is mostly dominated by two sources 

(GX 339—4 and V404 Cyg) observed in a single outburst. In this paper, we report 
on a series of radio and X-ray observations of the recurrent black hole GX 339—4 
with the Australia Telescope Compact Array, the Rossi X-ray Timing Explorer and 
the Swift satellites. With our new long term campaign, we now have a total of 88 
quasi-simultaneous radio and X-ray observations of GX 339—4 during its hard state, 
\Q ■ covering a total of seven outbursts over a 15-year period. Our new measurements rep- 

resent the largest sample for a stellar mass black hole, without any bias from distance 
uncertainties, over the largest flux variations and down to a level that could be close 
to quiescence, making GX 339—4 the reference source for comparison with other ac- 
1 creting sources (black holes, neutrons stars, white dwarfs and active galactic nuclei). 

Our results demonstrate a very strong and stable coupling between radio and X-ray 
emission, despite several outbursts of different nature and separated by a period of 
quiescence. The radio and X-ray luminosity correlation of the form Lx oc 01 
confirms the non-linear coupling between the jet and the inner accretion flow powers 
r jjj " and better defines the standard correlation track in the radio- X-ray diagram for stellar 

& \ mass black holes. We further note epochs of deviations from the fit that significantly 

exceed the measurement uncertainties, especially during the time of formation and de- 
struction of the self absorbed compact jets. The jet luminosity could appear brighter 
(up to a factor 2) during the decay compared to the rise for a given X-ray luminos- 
ity, possibly related to the compact jets. We furthermore connect the radio/X-ray 
measurements to the near- infrared/X-ray empirical correlation in GX 339—4, further 
demonstrating a coupled correlation between these three frequency ranges. The level 
of radio emission would then be tied to the near-infrared emission, possibly by the 
evolution of the broad-band properties of the jets. We further incorporated our new 
data of GX 339—4 in a more global study of black hole candidates strongly supporting 
a scale invariance in the jet-accretion coupling of accreting black holes, and confirms 
the existence of two populations of sources in the radio/X-ray diagram. 

Key words: accretion, accretion discs - binaries: general - ISM: jets and outflows - 
radio continuum: stars - X-rays: stars - stars: individual (GX 339-4). 
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Stellar mass black holes (BH) in accreting binaries un- 
dergo occasional outbursts with transitions between differ- 
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ent spectral states, usually defined by their X-ray spec- 
tral and tim ing properties (|McClintock fc Remillardl 120061 ; 
lBelloni|[2010h . In addition to the behaviour of the inner ac- 
cretion flow, outflows are also crucial for the chara cterisa- 
tion of these different spectral states (|Fenderll2006l ). espe- 
cially in light of the strong coupling between the jets and 
the disc. Two main spectral states can be defined as follows: 
1) the hard state is characterised by a comptonised X-ra 
spectrum with a weak or absent (e.g. iTomsick et al.ll20Q; 
thermal contribution from the accretion disc, and the pres- 
ence of powerful self-absorbed compact jet s JCorbel et alj 
200d; |Ph awan. Mirabel fc Rodriguez! l200d ; iFenderl l200lF 



Stirling et al .1120011 ). 2) the soft state is mostly dominated by 
the thermal emissio n from the accretion disc and the absence 
of relativistic jets (|Fender et al J 1 19991 ; ICoriat et aT1l2011bl ; 
I Russell et al.|[201ll ). We further note the intermediate states 
that are associated wit h many chan ges in the X-ray spectral 
and timin g properties ( Belloni 201 0]) and the inflow - outflow 



coupling (ICorbel et al 



.20041; iFender. Belloni fc Galldl2004 
iRodriguez et al.ll200Sl ; Fender. Homan fc Bellonill2009l V 

Weakly accreting BHs in the hard state have been the 
focus of many recent studies due to the presence of sev- 
eral emission processes that are very difficult to disentangle 
in this regime. Indeed, three main components are thought 
to dominate the energetic output of the system: the stan- 
dard optically thick and geometrically thin accretion disc 
(even tually truncated in its inner part) dShakura fc Sunvaevl 
1973h. the s elf-absorbed compact jets (|Corbel et al. 200d ; 



Fenderll200ll), and an optica lly thin corona of hot elec trons 
( Sunvaev fc Titarchuk|[l98d ; see also lXie fc Yuanl2012l ) that 



may or may not be the base of the compact jets (e.g. 
iMarkoff. Nowak fc Wilms! 120051 ). Emission from the com- 
panion star is usually negligible for low mass X-ray binaries 
in outburst. 

With the development of large campaign of multi- 
wavelength observations, several attempts have been 
conducted to model the spectral energy distribution 
(SED) of BHs in their regimes o f faint X-ray emis- 
sion (e.g. IMcClintock et all l200ll ; IMarkoff et all 120031 ; 
I Yuan. Cui fc Naravanl 1200 5l ). It turns out that optically 
thick synchrotron emission from the compact jets domi- 
nates the radio to near-infrared range in the hard state 
(ICorbel fc Fended200j;lHynes et al.l2003|;lBuxton fc Bailvn 
2004bl ; iHoman et al.ll2005l ; iRussell et al.H2006l ; ICoriat etal 



20091 ). In the hard X rays, inverse Comptonisation of the 
ambient photon field is the main emission process, with 
an active debate on whether it is thermal or non-thermal 
Com ptonisation, or synchrotron self-Compton (see f or a re- 
view IMcClintock fc Remillardl l200d ; iMarkofj l20ld ). Opti- 
cally thin synchrotron emission f rom the compact jets may 
also play a role at high energy (IMarko ff. Falc ke fc Fender 
2001 ; I Homan et al. I l2005l ; IRussell et al.ll20ld~ l Laurent et al 



20111 ) 



Another way to probe the emission properties of 
these BHs systems is to search for empirical relations 
between various wavelengths. The tight non-linear power 
law correlation between radio an d X-ray emission (of the 
form L Rad oc £x ' 7 in G X 339-4 (|Hannikainen et al.|[l99g| ; 
I Corbel et al1l200d . 120031 ) demonstrated the strong coupling 
between the compact jets and the X-ray emitting media. 
This correlation originally brought to light the possibility of 
optically thin synchrotron X-ray emission from the compact 



jets of GX 339-4 (|Corbel et al.ll2003l ; IMarkoff et aill2003l ) 
or alternatively from a radiatively inefficient a ccretion flow 
(iMerloni. Heinz fc di Matteoll2003l ; iHein j|2004 but see also 
IKording. Falcke fc Corbell |2006| ). iGallo. Fender fc Poolevl 
|2003l) extended this correlation study to a larger sample 
of BHs and proposed a universal correlation (still with an 
index of ~ 0.7) between the radio and X-ray luminosity for 
Galactic BHs. The correlation observed in the hard state 
appears to be maintained down to the quiescent level of 
at least two so urces (V404 Cyg, A 0620—00 and possibly 



also GX 339-4: ICorbel et al.ll2003l; IGallo. Fender fc Poolevl 



120031 ; IGallo et al.ll2006l ; ICorbel. Koerding fc Kaaretll2008l ). 

By taking into account the mass of the BHs, this correlation 
was extended to active galactic nuclei and led to the 
definition of the fundamental plane of BH activity and 
the possibility of univer sal scaling laws for accreting BHs 



Falcke. Kording & Markofl 


120041; IKording, Falcke & Corbell 


20061; Wang, Wu & Kong 


20061; Giiltekin et all 120091; 


Plotkin et al. 201ll). 



However, it is important to note that the 
universal radio/X-ray c orrelation presented by 
IGallo. Fender fc Poolevl (|2003h for BH binaries is dominated 
by two sources (GX 339-4 from ICorbel et all |2003l ) and 
V404 Cyg) with few measurements from additional sources 
that are often close to transition to the softer states. As the 
fundamental plane relies on the correlation for the Galactic 
BHs, it is important to verify how solid both the Galactic 
sampl e and the correlation are. I Corbel. Koerding fc KaaretJ 
(|2008l ) revisited the correlation for V404 Cyg and confirmed 
the tight correlation from outburst down to quiescence 
with an index of ~ 0.6 that could be consistent with X 
rays emanating from synchrotron self-compton emission 
at the base of the compact jets or Comptonisation from 
an inefficient accretion flow. A similarly tight correlation 
was also found between X-ray and optical/near- infrared 
(hereafter OIR) emission in a sample of BH candidates 
BHC) and also in GX 339—4 over different outburst s 
Homan et~aill2005l ; IRussell et al]|200d ; ICoriat et al.ll2009l ). 
Alternatively to jet or ADAF models, synchrotron emission 
from an hybrid electron distribution in a hot accretion flow 
has also been invoked to explain th e OIR/X-ray correlation 
jVeledina. Poutanen fc Vurmll201ll ). 

In recent years, significant efforts have been un- 
dertaken to observe most new BHs transients in out- 
burst at radio frequencies. This led to the discov- 
ery of several sources lying outside the scatter (hence 
their names "outliers") of the original radio/X-ra y cor- 
relation (e.g. XTE J1650-500. ICorbel et all |2004| ; IGR 



J17497-2821. IRodriguez et al.l [20071 Swift J1753.5-0127, 



ICadolle Bel et all 120071 ; ISoleri et all l20ld ). These out 
liers lie significantly below the standard correlation 
jGallo. Miller fc Fender! |2012[ ). Usin g a large sa m ple for 
the BH transient H 1743-322, ICoriat et all |2011bl ) 
found a steeper correlation index of ~ 1.4 for the out- 
liers (and possibly also other sources l ike Cyg X— 1 
(IZdziarski et al.l 1201 ll ) and GRS 1915+105 (|Rushton et al.l 
2010h). At fainter fl u xes, these out l iers (at least H 1743-322, 



Jonker et al.l (|20ld ) ; ICoriat et all |2011bl ) and p ossibly now 



also MAXI J 1659- 152 and XTE J1752-223: iRatti et all 
|2012| ; Ijonker et al.l |2012| ) may return to the standard cor- 
relation. It remains unclear how these outliers are related to 
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typical black holes such as GX 339-4 or V404 Cyg. With a 
lobah study of H 1743—322 and other sources. ICoriat et al.l 
2011bl ) proposed that these outliers may represent a large 
subset of BHs radiating efficiently in the hard state. This 
is a change of paradigm for BHs in the hard state, usually 
believed to be in a regime of inefficient accretion (as also 
implied by the standard radio/X-ray correlation). Alterna- 
tively, a di fferent coupling betw een jet power and mass accre- 
tion rate (|Coriat et al.ll2011bl) may also explain the steeper 
correlation index (see also Soleri fc Fenderll201ll ). 



In this paper, we review the properties of the correla- 
tion for GX 339—4 based on a very large sample of new 
data acquired during several outbursts in the last decade. 
Indeed, amongst the Galactic sources in the standard cor- 
relation (in fact dominated by one outburst for GX 339—4 
and V404 Cyg), none has been observe d over several out - 
bursts. The original work on GX 339-4 (|Corbel et al.ll2003l ) 
included data from the 1997 extend ed hard state (e.g. 
lHarmon et al. 1 11994 ICorbel et al1l2000l ) and the decay from 
the 1998-99 outburst. The inclusion of two radio observa- 
tions from the rise of the 2002-03 outburst pointed out that 
GX 339—4 could fo llow another track for different outbursts 
jNowak et al.ll2005h . 

Due to its recurrent outburst activity, GX 339—4 is an 
ideal target to probe the evolution of the correlation over 
distinct outbursts. After the 1998-99 outburs t, GX 339-4 
was found in quiescence in 2000 and 2001 (|Corbel et alj 
120031 ). Since then, seven outbursts were observed: 2002- 
03 dHoman et al.l 120051 ; iBelloni et ail 120051). 2004-0 5 
jBelloni et al.ll2006l ), 2006, 2007 (e g. iTomsick et al.ll2008l) 
2008 2009, 2010-11 (IShidatsu et al.ll201ll;ICadolle Bel et al l 
l201ll ; iMotta et al.| [2011: Rah oui et al1l2012l) and a tempo- 
rary reactivation in 2012 llLewis Russell fc Shahbad lioil ; 
iMaccarone. Russell fc Lewis! I2012T ). GX 339-4 is believed 
to harbour a BH with a ma ss > 5.8 Mp, in a system with 
a low mass star companion (|Hvnes et al. l2004h. GX 339— 4 



is loc ated at a favoured distance of 8 kpc (Zdziarski ct al. 
l2004h . In section 2, we present the radio and X-ray observa- 
tions used in this analysis. The results for selected outbursts, 
as well as the re-evaluation of the radio/X-ray correlation 
using all data, are presented in Section 3. We then discuss 
these results in Section 4 in light of the multi-wavelength 
behaviour of GX 339—4. Our conclusions are summarised in 
Section 5. 



2 OBSERVATIONS AND DATA ANALYSIS 
2.1 ATCA observations 

All radio observations discussed in this paper were con- 
ducted with the Australia Telescope Compact Array 
(ATCA) located in Narrabri, New South Wales, Australia. 
The ATCA synthesis telescope is an east- west array consist- 
ing of six 22 m antennas. It uses orthogonal linearly polarised 
feeds and records full Stokes parameters. We carried out all 
observations simultaneously at 4.8 GHz (6.3 cm) / 8.64 GHz 
(3.5 cm) and in a limited cases also at 1.384 GHz (21.7 cm) / 
2.496 GHz (12.0 cm). The ATCA was upgraded in April 2009 
with th e new Compact Arr ay Broadband Backend (CABB) 
system (|Wilson et al.ll201lT ) resulting in significant improve- 
ment in sensitivity (an increase in bandwidth from 128 MHz 



to 2 GHz). This implied a small change in central frequen- 
cies and the most recent observations were carried out at 
the frequencies of 5.5 and 9 GHz. Various array configura- 
tions have been used during these campaigns. All observa- 
tions have been conducted by us, with the exception of four 
observati ons (PI: M. Rupe n) in March 2004 (briefly men- 
tioned in lMiller et al.|[2006T l. To compare with the previous 
outbursts, we also u sed the orig inal data of GX 339 —4 from 
ICorbel et all (|2005t ) updated in lNowak et~aTI (120051 ). 

The amplitude and band-pass calibrator was 
PKS 1934—638, and the antennas' gain and phase 
calibration, as well as the polarisation leakage, were usually 
derived from regular observations of the nearby (less than 
a degree away) calibrator PMN 1646—50. The editing, 
calibration, Fourier transformation with multifrequency 
algorithms, deconvolution, and image an alysis were per- 
form ed using the MIRIAD software package (|Sault fc Killeenl 
ll99Sft . In short observations with bad coverage of the u-v 
plane, flux densities were measured directly in the u-v 
plane with MIRIAD task uvfit and consequently checked 
afterwards with direct imaging of simulated datasets. We 
only present in this paper the data when GX 339—4 was 
in the hard state according to the X-ray behaviour (black 
squares in the vertical track of the Hardness-Intensity 
Diagrams (Fig. [1} . 



2.2 RXTE observations 

The Rossi X-ray Timing Explorer (RXTE) has conducted 
(almost daily) monitoring observations of GX 339—4 dur- 
ing almost a l l out b ursts discussed here (see for example 
Homan et ail 120051: IBelloni et al.l IgOOfj; IMotta et all |201ll : 



Stiele et al 1l201ll : ICadolle Bel et al.ll20nl ), with an exposure 
time of typically 1 to 3 ks for each observations. 

We performed spectral analysis using data from the 
Proportional Counter Array (PCA, Ijahoda et"ail 120061 ) 
and the High Energy X-ray Timing Experiment (HEXTE, 
iRothschild et all 1 19981 ) from observations when GX 339-4 
was in the hard state. The data were reduced using HEASOFT 
software package v6.11 following the standard steps de- 
scribed in the (RXTE) cookboolo We extracted PCA spec- 
tra from the top layer of the Proportional Counter Unit 
(PCU) 2 which is the best calibrated detector out of the 
five PCUs and the only one operational across all the obser- 
vations we analysed. We produced the associated response 
matrix and added a systematic uncertainty of 0.5 % to all 
spectral channels to account for PCA calibration uncertain- 
ties. In addition, we used the background model (faint or 
bright) appropriate to the brightness level of GX 339—4 to 
create background spectra. 

For HEXTE, we produced a response matrix and ap- 
plied the necessary dead-time correction. The HEXTE back- 
ground is measured throughout the observation by alternat- 
ing between the source and background fields every 32 s. The 
data from the background regions were then merged. When 
possible we used data from both detector A and B to extract 
source and background spectra. However, from December 
2005, due to problems in the rocking motion of Cluster A, 
we extracted spectra from Cluster B only. HEXTE channels 



1 http:/ /heasarc. gsfc.nasa.gov/docs/xtc/data_analysis. html 



4 S. Corbel et al. 




Hardness Ratio 



GX 339—4 was detected at its faintest level of radio emis- 
sion. 

In addition, we constructed Hardness Intensity Dia- 
grams (HIDs, see Fig.[T|) from PCA data. We produced back- 
ground subtracted lightcurves binned at 16 s using PCA 
Standard 2 mode data from the PCU2 (all layers). The 
lightcurves were divided into three energy bands, 2.5-6.1, 

9.4- 18.5 and 2.5-18.5 keV (according to the correspondence 
between fixed energy channels and energy valid for the PCA 
gain epoch 5). We defined the hardness ratio (HR) as the 
ratio of count rate in the bands 9.4-18.5 keV and 2.5-6.1 
keV and the intensity as the count rate in the energy band 

2.5- 18.5 keV. 



Figure 1. B.XTE/PCA hardness intensity diagrams (HID) cov- 
ering all outbursts of GX 339-4 in the 2000-2011 period (see end 
of section 12.21 for definition of the X-ray hardness and intensity) . 
The squares highlight the epoch of our radio observations (not 
including the ones conducted at very low flux) . Those on the ver- 
tical track are observations purely in the hard state, whereas the 
few ones on the horizontal branches correspond to period of com- 
pact jet destruction in 2010, and formation in 2007 and 2011 (see 
Section. [374t. 



were grouped by 4 due to the low count rate in most of the 
observations. On Dec. 14 th 2009, Cluster B stopped rocking 
as well. From this date, we thus used only PCA data in our 
analysis. 

The X-ray fluxes in various energy bands were de- 
termined by fitting simultaneously the PCA and HEXTE 
spectra in XSPEC V12.7. We used a model consist- 
ing of a multi-color disc blackbody(ezdisA;66, not required 
at low flux), a power-law with a high energy cut-off 
when required (power x highecut), and a Gaussian emis- 
sion line (gaussian) constrained in energy between 6 
and 7 keV. Interstellar absorption was taken into ac- 
count using the phabs model w i th cr oss-sections from 
Balucinsk a-Church fc McCammonl (Il992f) and abundances 
from I Wilms. Allen, fc McCravl |2000r i. We fixed the hydro- 



gen column density to Nu = 6 x 10 21 c m" 2 , which is a 
typical value observed in G X 339—4 (e.g. IZdziarski et alj 
l2004l ; [c"adolle Bel et al.l201ll ). A multiplicative constant was 
added to the model to allow for differences in the over- 
all normalisation between both instruments. All fluxes pre- 
sented have been normalised to PCA and have been cor- 
rected for the interstellar absorption (its effect is almost 
negligible for GX 339—4 above 3 keV). The flux error bars 
were chosen to be the largest of the statistical error, or 
3%, which is a good estimate of the PCA systematic un- 
certainty. Furt hermore, the X-ray fluxes from the origi- 



nal correlation (|Corbel et al.ll200a ) have also been updated 
with the last version of the RXTE soft ware and calibra- 
tion. The level of Galactic Ridge emission jRevnivtsev et al.l 
l200d : iRevnivtsev fc Sazonovll2007l ; lEbisawa et al.ll200St ) has 
been estimated from long PCA observations taken when 
GX 339—4 was in quiescence (and th erefore below PCA de- 
tection level, see ICoriat et all 12009] for details). This level 
of 4 x 10~ 12 erg s -1 cm -2 (3-9 keV band) has been sub- 
tracted from all X-ray fluxes. Quantifying this level was im- 
portant for the discussion on the level of radio emission when 



2.3 Swift observations 

We used measurements made by the Swift XRT instrument 
to extend this study to low flux levels. Our analysis included 
one 3.0 ks observatio n made on 2009 Oc tober 30 that was 
previously reported bv lYen fc Konel (|2009T l , five observations 
made in 2011 on April 13 (1.2 ks), April 27 (4.9 ks), May 
11 (4.8 ks), and September 9 (2.1 ks), and two observations 
made in 2012 on June 23 (0.15 ks) and June 26 (0.43 ks). 
Although more Swift observations are available near these 
times, these are the observations that are closest in time to 
our radio observations. 

For all eight of these observations, XRT was operated in 
photon counting (PC) mode, providing 2-dimensional imag- 
ing. We analyzed the data using the HEASOFT software 
package. For the 2009 and 2011 observations, we reprocessed 
the data to produce cleaned event lists using the routine 
xrtpipeline, and for the 2012 observations, we used the 
output of the standard pipeline processing. We used xselect 
to extract source spectra from a circular aperture with a 20 
pixel (47") radius and background spectra from an annular 
region centered on the source. After background subtraction, 
the 0.5-10 keV count rates for the five observations range 
from 0.007 ± 0.002 s _1 to 0.14 ± 0.02 s~\ which are low 
enough for photon pile-up to be negligible. We used the PC- 
mode response file swxpc0tol2s6_20010101v012.rmf, and the 
routine xrtmkarf to produce the ancillary response ("arf") 
files. For the spectral fits described below, we combined the 
2012 June 23 and June 26 spectra due to the short exposures 
and the fact that the radio observation was made on June 
24th in between the two Swift observations. 

To determine the X-ray flux, we used the XSPEC soft- 
ware package to fit the energy spectra with a power-law 
model with interstellar absorption. We used Cash statistics 
for the fits, which does not require binning of the spectra. 
We fixed the column density to Nh = 6 x 10 21 cm -2 , and 
the photon indexes obtained for the non-quiescent observa- 
tions ranged from T = 1.5 ±0.3 to 2.0 ±0.3 (90% confidence 
errors) . For the three quiescent observations, the photon in- 
dex was not well-constrained, and we fixed it to T = 2. The 
main goal of these spectral fits is to obtain the unabsorbed 
flux for each spectrum, and we accomplished this by using 
the model pegpwrlw, which has flux in an adjustable energy 
band as a free parameter, to represent the power-law. 
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Figure 2. RXTE/ASM soft X-ray (1.5-12 keV) light-curve of 
GX 339-4 for all RXTE lifetime (January 1996 to December 
2011). The major outbursts, that include a full transition to 
the soft state, are indicated by their epoch. In addition, the ar- 
rows point to additional fainter outbursts that are not visible 
with ASM (but see the Swift light-curve on Figure [3} because 
GX 339—4 was only in the hard state. 



2.4 Notes on some individual observations 

Almost all radio and X-ray observations were conducted 
quasi simultaneously (less than a few hours). However, in 
only very limited cases, we interpolated (usually with a sec- 
ond degree polynomial on a linear scale) the monotonously 
rising/decaying X-ray flux in order to get the most precise 
flux at the time of the radio observations. 

One single ATCA observation was conducted in 2006 
on November 5 (MJD 54044.2) and it occurred well after 
the end of the faint 2006 outburst (from January to June 
2006, see the Swift/BAT ligthcurve in Fig- [3j . Inspection of 
the RXTE/PCA light-curve (not shown) revealed that the 
reported 2007 outburst already started, at least, as early as 
2006 November 12 (date of the first available PCA obser- 
vation) with an 3-9 keV unabsorbed X-ray flux of 5.56 x 
10 -11 erg s _1 cnT 2 . Thirteen RXTE/PCA observations oc- 
curred in a thirteen-day period and it showed an X -ray flux 
that was steadily rising (see also lSwank et al]|2006T ). A fit to 
the rising X-ray flux leads to an extrapolated 3-9 keV flux 
of 2.2 x 10 -11 erg s _1 cm~ 2 (with a very conservative error 
of 30%) for 2006, November 5 (i.e 5 days earlier). 

During the 2008-2009 period, GX 339-4 is mostly de- 
tected by the hard X-ray all sky monitor (aka Swift/SAT), 
indicative of hard state activity only. For the observation on 
MJD 55002, the X-ray flux was taken from the RXTE/PCA 
Galactic Bulge Monitoring of Craig MarkwardlQ. On MJD 
54469, the X-ray flux was estimated from an average of the 
RXTE/PCA observations conducted 6 days before and af- 
ter, at a period when GX 339—4 was not very variable. Fi- 
nally, on 2011 May 11, we used the X-ray flux derived from a 
Chandra observations (to be discussed elsewhere) conducted 
at low flux. 

Regarding the radio observations, we note that in a few 
cases, it was also necessary to take into account the presence 



http: / /asd. gsfc.nasa.gov/Craig.Markwardt/galscan/main. html 



Figure 3. Swift/BAT hard X-ray (15-50 keV) light-curve of 
GX 339—4 corresponding to Fig. [2] after Swift launch (February 
2005 to December 2011). Outbursts are identified by their epoch. 
We note that after the 2007 outburst, GX 339—4 brightened again 
(noted '2007 reflare') for few months and stayed only in the hard 
state, as in 2006, 2008 and 2009. 



of an additional component associated with the interaction 
of pr eviously ejected material with the inters tellar medium 
(e.g. iGallo et all |2004| ; ICorbel et all l2010bl ) similarly to 
what has been detected in XTE J1550-564 (|Corbel et all 
|2002[ ). However, this did not concern any hard state mea- 
surements, but only the observations during the re-ignition 
of the compact jets in the 2007 soft to hard state transition. 
Finally, the first ATCA observations in 2002 was only made 
at 1.384 and 2.496 GHz and flux densiti es have been extrap - 
olated to 8.6 GHz assuming a typical l|Corbel et al.l I2003T ) 
radio spectral index a =+0.12, taking a definition of the 
flux density, S„, as S v cx v a (consistent with the 1.384 to 
2.496 GHz spectral index of 0.11 ± 0.11). Similarly, on 2002 
April 5, the flux at 8.6 GHz is deduced from the detection 
at 4.8 GHz (no data usable at 8.6 GHz for this observa- 
tions). All our measurements representing the whole sample 
of GX 339—4 in the hard state are tabulated in Table 1. 



3 RESULTS 

3.1 A long series of outbursts in a decade 

In the last decade, GX 339—4 has been one of the most ac- 
tive accreting BHs with many outbursts (Figs. [ 2] and O. The 
work on the radio/X-ray flux correlation by Corbel et alj 
(|2003h used data from 1997, when GX 339—4 was in a period 
of extended and persistent hard state, and data from the de- 
cay of the 1998/99 outburst. Four major outbursts have been 
observed since 2000: 2002/03, 2004/05, 2007 and 2010/11. In 
addition, the analysis of the Swift/BAT light-curve reveals 
weaker outbursts in 2006, 2008 and 2009, when GX 339-4 
was observed only in the hard state. After the decay of 
the 2007 outburst, GX 339—4 brightened again (denned as 
"2007 reflare" in Fig. [3j but stayed only in the hard state. 
We have plotted in Fig. [T] the full hardness intensity dia- 
gram (HID). For 1998/99, very few RXTE observations were 
conducted and therefore these HIDs are not represented. 
For a general and recent overview of the activity, see e.g. 
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Table 1. Summary of data used for the radio/X-ray flux correlation in GX 339—4 



Calendar Date 


MJD 


Flux density 






Radio 


X-ray flux (3-9 keV) 


X-ray 


Comment 


CUT) 




4.8 or 5.5 GHz 


8.6 or 


9.0 GHz 


spectral 


(in unit of 10" 10 


hardness 








(mJy) 






index 


erg s — 1 cm -2 ) 






1997 Feb. 4 


50483.04 


N.A. 


9.10 


± 


0.10 


N.A. 


9.66 ± 0.34 


0.71 


C03, N05 


1997 Feb. 10 


50489.79 


N.A. 


8.20 


± 


0.20 


N.A. 


8.48 ± 0.30 


0.72 


C03, N05 


iyy i rcu. 1/ 


t^O/IOfi CI 

DU4yo.oi 


N.A. 


8.70 


± 


0.20 


TS.T A 


oil -t n 90 
o.ii it u.zy 


O 79 
U. / Z 


L/Uo, 1NU0 


iyyy rcu. iz 


K~\ 991 Ofi 

oizzi .yo 


6.34 ± 0.08 


4.60 


± 


0.08 


o kz: _i_ n r\A 

— U.00 it U.U4 


a qo -i- n i a. 
4.oy it u.io 


O K*7 
U.O / 


uUo, 1NUO, INa 


iyyy iviar. o 


til 9/10 71 
01Z4U. i 1 


6.07 ± 0.06 


5.74 


± 


0.06 


O 1 O -1- O OQ 
— U.1U it U.Uo 


A Qc; _|_ n i c 
4.o0 it U.IO 


O 1« 
U. / 


OUo, 1NU0, INa 


1999 Apr. 2 


51270.56 


4 7k -L n f|6 

rl. 1 O _1_ U.UU 


5.10 


± 


0.06 


+0.12 ± 0.04 


4.42 ± 0.16 


0.53 


C03, N05 


1999 Apr. 22 


51290.58 


2 qo + 06 

— ■ ■ .' — _i_ w . Kjyj 


3.20 


± 


0.06 


+0.16 ± 0.05 


2.05 ± 0.07 


0.55 


C03, N05 


1999 May. 14 


51312.65 


i 25 + 05 


1.44 


± 


0.06 


+0.24 ± 0.10 


0.72 ± 0.03 


0.52 


C03, N05 


1999 Jun. 25 


51354.69 


W. J-Tl _1_ U.OO 


0.24 


± 


0.05 


+0.92 ± 0.51 


(3.93 ± 0.28) x lO- 2 


N.A 


C03 


1 OOO Tul 7 

iyyy jui. ( 


^1 Qfifi ^fi 

OloOO.OO 


< 0.13 


0.12 


± 


0.04 


IN .A. 


f*i -i- n 1 o\ v 1 n - 2 
l^l.Oo it u.iyjxiu 


KT A 
IN . I\ 


L/Uo 


1 OOO Anir 17 

iyyy ivug. i 1 


c;i A 07 Q X, 
D14U / .oO 


N.A. 


0.27 


± 


0.07 


TS.T A 


{ o n -t n '^^^/1n— 2 

(^Z.U it U.o ) X 1U 


M A 
IN ./V 


t-Uo 


1 OOO Car. 1 

iyyy ocp. i 


K~\ /1 99 
014ZZ.OO 


N.A. 


0.32 


± 


0.05 


NT A 


/ 1 ok -I- n inWin - 2 
^i.yo it u.iujxiu 


NT A 
IN .I\ 


t-Uo 


1 OOO Q/-ii-x i 

iyyy ocp. 10 


- , ^1 A Qfi 
~ D14oO 


N.A. 


< 


0.06 


NT A 
IN 


/'k 7 _i_ 1 n\vin-3 

\p. ( it l.UJ X 1U 


TJ A 
IN .I\ 


L/Uo 


zuuz /\pr. 4 


f^9QfiO CO 

ozooo.uy 


N.A. 


5.95 


± 0.15 


i n 11 _|_ n 1 1 
4-U.ll it U. 11 


1 fi9 4- H (i1 

lo.OZ it U.Ol 


O KQ 

U.Oo 


1NU0, Q|Z. 1| 


2002 Apr. 5 


52369.66 


fi 77 + 1 5 

U . 1 1 _1_ VJ . J. ij 


7.26 


± 0.20 


N.A. 


22.86 ± 0.73 


0.57 


S\2A\ 


2002 Apr. 7 


52371.70 


8 39 + 06 


8.27 


± 0.07 


-0.02 ± 0.02 


30.63 ± 0.96 


0.54 




9009 A i-ii- 1 o 
zuuz /ipr. lo 


COQQO Of? 


i q i q 4. n 


14.27 ± 0.05 


i n 1 q _i_ n m 
H-U.lo it U.U1 


c:q Ofi -U 1 
Oo.oO it 1.00 


o c;i 
U.Ol 


1NU0 


OflflQ T\/Toi r Oc; 

zuuo iviay zo 


C07Q/1 ni 

DZ / S4.U1 


n fjn 4. n n^i 

U.UU _1_ W.L/iJ 


0.77 


± 


0.06 


I n A 9 -l- O 1 O 

-|-U.4Z it u.iy 


n i Qfi _i_ n mn 
U.loO it U.UoU 


n fin 
U.OU 




900/1 T?oK /I 

ZUU4 reD. 4 


c^QOQO 07 


N.A. 


0.70 


± 


0.04 


T\I A 
IN .A. 


n ik -4- n n 1 !! 

U.140 it U.Uol 


n 7n 
U. / U 




900/1 T7V>K 1 Q 

ZUU4 reu. lo 


OOU4o.Oy 


1.27 ± 0.07 


1.13 


± 


0.08 


90 4- n 1 ^ 
— u.zu m u.io 


n 7/1Q J_ 097 

u. / 4y it u.uz / 


^9 

u. oz 




Oflfl/I 17nK 9/1 
ZUU4 rGD. Z4 


c;qo^o no 
oouoy.uy 


1.77 i 0.15 


1.84 


± 


0.20 


i n ri7 _i_ n 9Q 
-|-U.U/ it U.Zo 


9 Qfi -U n oo 
Z.oO it U.Uo 


n fi9 

U.OZ 




900A A/Tar- 1 fi 
ZUU4 IViar. ID 


oouou. i y 


4 7K 4. n ri6 

rl. 1 _1_ U.UU 


4.88 


± 


0.06 


in 10-I-0 07 
-T"U. 1U III U.U ( 


Q 07 -1- O Q9 

y.ui iii u.oz 


fil 
U.Ol 




2004 Mar. 17 


53081.79 


4.92 ± 0.05 


4.84 


± 


0.11 


-0.03 ± 0.05 


9.55 ± 0.34 


0.61 




2004 Mar. 18 


53082.79 


5.15 ± 0.05 


4.98 


± 


0.11 


-0.05 ± 0.05 


9.92 ± 0.35 


0.61 




onn/i A/T-iv i o 
ZUU4 iviar. iy 


[ronQQ 70 

Ov5Uoo. / y 


5 17-1-0 04 


5.20 


± 


0.10 


i n m _L n r\A 
+U.U1 it U.U4 


1 n nfi _l n q c; 
1U.UO it U.oo 


n fin 
U.OU 




2005 Apr. 21 


53481.73 


4 qk 4. n 05 


4.73 


± 


0.05 


+0.14 ± 0.03 


3.75 ± 0.13 


0.55 




2005 Apr. 24 


53484.67 


4 39 + 06 


4.23 


± 


0.08 


-0.06 ± 0.04 


2.77 ± 0.10 


0.57 




OOO^ A iAi- OC 

ZUUO Apr. Zo 


^Q/IOO t^Q 
Do4oo.0o 


q n7 4- n 1 n 


3.46 


± 


0.13 


n on -U o no 
Tu. ZU it U.Uo 


i oc; _i_ n 07 
l.oO it U.U/ 


n c;o 
U.Oo 




900^ A i-ii- 90 

zuuo Apr. zy 


c;'?/! oo fit; 
Do4oy.00 


2 ss 4- n 08 


3.32 


± 


0.10 


i n O/i _i_ n fi7 
-|-U.Z4 it U.U/ 


i -U n ofi 
1.00 it U.UO 


n c;o 
U.Oo 




onn^ A i-ii- qo 
zuuo -rt-pr. ou 


oo4yu.Uo 


2 c:'? 4- n 05 


2.94 


± 


0.07 


-tU.ZD it U.UO 


1 -U O Ofi 
1.0/ it U.UO 


n c;o 
U.Oo 




OOO^ l\/fr,,, O 

zuuo iviay o 


0o4yo.oo 


i 7Q -L n in 

-L . ( O ^ U. J-U 


1.92 


± 


0.14 


j_ n 1 o -U n i a 
-f-U-lo it U.lu 


1 17 i n o/f 

1.1 / it U.U4 


n c;o 
u.oy 




900^ l\/Toir A 

zuuo iviay 4 


^Q/IO/I cc 
OD4y4. oo 


i Q7 _|_ n no, 

l.ul _1_ W.L/O 


1.99 


± 


0.10 


i n 11 _i_ n 1 1 
-(-U.11 it U.ll 


i 9^ -U n oc; 

1.Z0 it U.UO 


n c;o 
U.Oo 




900^ AAoi, fi 

zuuo iviay o 


^Q/IOfi 00 


1.42 ± 0.09 


1.69 


± 


0.12 


i n on -l- n i f\ 
tU. oU it U.lu 


n ofi -U n o/i 
u.yo it U.U4 


n c;o 
U.Oo 




900^ A/Toxr 1 9 

zuuo iviay iz 


OoOUZ.Qo 


1.08 ± 0.12 


1.00 


± 


0.18 


n 1 q -U n Qfi 
— U.lo it U.oU 


n 7Q -4- n HQ 

U. (O ± U.Uo 


n c^fi 
U.OO 




900fi 1ST™ 7 ^ 
ZUUD 1NOV. D 


c;/i O/i /i i 7 

04U44. 1 ( 


w.iji _i_ w . yj\j 


0.60 


± 


0.07 


i n 09 -i- n ia 
-+-U.OZ it U.o4 


n 99 -U n oo 

U.ZZ it U.Uo 


INT A 
IN .A. 


c\n A\ 


2007 Feb. 4a 


54135.10 


19 5 + 03 


17.0 


± 1.0 


-0.23 ± 0.10 


54.84 ± 1.92 


0.50 


S|3.4.2| 


2007 Feb. 4b 


54135.77 


21.9 ± 0.3 


22.5 ± 0.3 


+0.05 ± 0.03 


56.80 ± 1.99 


0.50 


S|3.4.2| 


9007 A/Toir Q1 

zuu i iviay oi 


c;/19^1 7/1 
04Z01. / 4 


4.34 ± 0.14 


4.37 


± 


0.09 


in m -U o nfi 
-|-U.U1 it U.UO 


9 1 o -U n oo 

Z.1U it U.Uo 


n ^.i 

U.O / 




9007 Tun fi 

ZUU / JUII. D 


t^/1 9^7 7/1 
04Z0 / . / 4 


2.63 ± 0.15 


2.63 


± 


0.18 


j_ n on j- n i a 
-|-U.UU it U.14 


i ki -U n ok 

l.Oo it U.UO 


n c;o 
u.oy 




9007 Tnn 1 1 

zuu i j un. 1 1 


^/1 9fi9 7Q 
04ZDZ. f O 


1.77 ± 0.18 


2.01 


± 


0.15 


j_ n 99 _u n 91 
-f-u.zz it U.Zl 


i 90 -U n o^ 
l.Zo it U.UO 


n c;o 
U.Oo 




9007 Tnn OC; 

zuu i j un. zo 


X.A 97fi 71 
04Z ( D. ( 1 


1.69 ± 0.06 


1.69 


± 


0.05 


in nn -U n no 
-f-U-UU it U.Uo 


i 9fi 4. n nt; 
l.ZO it U.UO 


n c;o 
u.oy 




9007 Tim 90 

zuu / jun. zy 


t^/1 90O /I fi 
D4ZSU.40 


1.7 ± 0.2 


2.0 


± 


0.2 


i n oo -U n 9fi 
Tu. Zo it U.ZO 


i _i_ n n^ 
l.oo it U.UO 


n fin 
U.OU 




2007 Jul. 4 


54285.52 


1.7 ± 0.2 


2.1 


± 


0.2 


+0.36 ± 0.26 


1.81 ± 0.07 


0.61 




2007 Jul. 13 


54294.50 


2.36 ± 0.05 


2.66 


± 


0.05 


+0.20 ± 0.05 


2.50 ± 0.09 


0.62 




2007 Aug. 22 


54334.50 


3.30 ± 0.05 


2.95 


± 


0.07 


-0.19 ± 0.05 


3.67 ± 0.13 


0.62 




2007 Nov. 3 


54407.25 


0.55 ± 0.06 


0.80 


± 


0.07 


+0.63 ± 0.24 


0.25 ± 0.04 


0.57 




2007 Nov. 27 


54431.99 


0.59 ± 0.07 


0.48 


± 


0.07 


-0.35 ± 0.32 


0.25 ± 0.06 


0.52 




2008 Jan. 4 


54469.05 


0.70 ± 0.05 


0.65 


± 


0.05 


-0.13 ± 0.18 


0.275 ± 0.063 


N.A. 


sua 


2008 Jun. 26 


54643.77 


0.77 ± 0.08 


1.16 


± 


0.10 


+0.70 ± 0.23 


0.714 ± 0.075 


0.62 




2008 Jul. 5 


54652.70 


1.08 ± 0.06 


1.24 


± 


0.07 


+0.24 ± 0.13 


0.875 ± 0.092 


0.62 




2008 Jul. 16 


54663.70 


1.21 ± 0.07 


1.51 


± 


0.06 


+0.38 ± 0.12 


1.17 ± 0.12 


0.64 




2008 Aug. 18 


54696.60 


1.10 ± 0.10 


1.18 


± 


0.10 


+0.12 ± 0.21 


0.858 ± 0.090 


0.62 




2008 Sep. 29 


54738.02 


0.69 ± 0.10 


0.91 


± 


0.10 


+0.47 ± 0.31 


0.454 ± 0.49 


0.57 




2008 Oct. 10 


54749.37 


0.86 ± 0.07 


0.73 


± 


0.10 


-0.28 ± 0.27 


0.166 ± 0.021 


0.76 
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Table 1 - continued Summary of data used for the radio/X-ray flux correlation in GX 339—4 



Calendar Date 


MJD 


Flux density 


Radio 


X-ray flux ( 3-9 keV) 


X-ray 


Comment 


CUT) 




4.8 or 5.5 GHz 


8.6 or 9.0 GHz 


spectral 


(in unit of 10~ 10 


hardness 








(mJy) 


index 


erg s _1 cm -2 ) 






onno Lin on 
zuuy jun. zu 


Kcnno £.0. 
oouuz.oy 


9 92 + 03 


2 7K 4- n (13 


_i_n A a 4- n c\a 


ni 4- n t^i 


n fi^ 
U.OD 


rt.A. 1 Hi \jDLVL 


zuuy oct. ou 


C K 1 OA QQ 

00104. 00 


nnmninpn 1 1 TiT^ifir 

v_Jl_J!IlLJllll_'V_i. U 1/ 1) \_> 1 




KT A 


/o -J- n aWin - 3 
yz.z ± v.ojx lu 




own! / Arvi 


om n Ton 01 
zuiu Jan. zi 


x. t^OI *7 OO 

oozi / .yz 


5.08 ± 0.04 


5.05 ± 0.05 


n ni 4. n ni 
— U.Ul ± U.Uo 


fi no 4- n 01 
D.UZ It u.zl 


n fifl 
U.00 




zuiu reD. lo 


c^o/in m 

00Z4U.U1 


U.l 1 _1_ U.UU 


k qn + n 1 


n no 4- n n/i 
— U.uy ± U.U4 


Q OB 4- n OO 
0.Z0 it u.zy 


n fifi 
U.00 






oozoo.oy 


7 22 + 10 


7 30 4- 10 


in no 4- n c\a 

TU.UZ ZLL U.U4 


1 a nfi 4- n 4Q 
14. uo it u.^y 


u.00 




2010 Mar. 6 


55261.89 


q f)2 + 1 


q «n 4- n 05 


+0.13 ± 0.03 


16.41 ± 0.58 


0.62 




2010 Mar. 7 


55262.91 


a 24 + 05 


a fir 4. n 1 n 


-0.05 ± 0.03 


16.49 ± 0.58 


0.62 




om n A/Tcr 1 a 
zuiu iviar. 14 


oozoy.ou 


1018 + 010 

HMO _1_ U.1U 


11 qo + 1 


n 00 -U n hq 
-+-U.ZZ it U.Uo 


on ni 4. n 7n 

ZU. Ul it U. (U 


n fin 
U.DU 




om n A/Tcr 1 a 
zuiu iviar. id 


OOZ ( l.DU 


1 n 85 + 1 

i-\J.(J<J _1_ U.1U 


12 ni + n m 


1 n 01 -U n nt 
-|~U.Z1 it U.Uo 


01 01 -4- n 11 
zi.yi ± u.n 


n £.0. 
u.oy 




Om D A/fnr 01 

zuiu iviar. zi 


K£9'7K QO 
OOZ ( O.oZ 


1 3 76 + 10 

J_v>. 1 «J _1_ U. 1U 


1 5 45 + 06 


1 n oa -i- n no 
-|-U.Z4 it U.UZ 


OQ 1R 4. n 00 

zo. id it u.yy 


n £.7 
U.O / 




Om n l\(Tnr OA 

ZUIU IVlcLr. Z4 


ooz / y. / y 


15 56 + 05 


1 a e:q 4. n ne; 


1 n Qfi _i_ n m 
-|-U.oD it U.Ul 


q no -4- 1 1 a 
00. uy it lid 


n k; k: 

u.00 




om n i\/Tnr oa 
zuiu iviar. zo 


K.KOQQ 7Q 
OOZoo. ( 


1 Q 4-8 -1-0 10 

J-C/.rtO _1_ U.1U 


21 88 + n in 


1 n oa -L n no 
+ U.Z4 it U.UZ 


qo no -U 1 QQ 
00. UZ it l.oo 


u.00 




om n i\ t 1- qi 
zuiu iviar. oi 


OOZoO. I O 


22 68 + 05 


2^; qA 4- n n^i 


1 n ov _i_ n m 
+U.Z/ it U.Ul 


A n f^V -4- 1 A 
4U.0/ it 1.4Z 


u.00 




om n A tai- 
zuiu Apr. z 


oozoo.y 1 


21 qc; + n 05 




1 n oq _i_ n m 
~hU. Zo it U.Ul 


A Q no -4-1 t^1 

4o.uy it 1.01 


n f^Q 
U.Oo 




om n A 1-11- q 
zuiu /Apr. 


oozsy.y 1 


18 84 + 10 


21.11 ± 0.15 


1 n 0*3 -U n no 

■fl). Zo it U.UZ 


AA Qfl -U 1 t^c; 
44. oD it 1.00 


n £.q 
U.Oo 




Om fl A nr A 

ZUIU JA.pr. 4 


KCOQfl CO 

oozyu. oz 


21 60 + 04 


23 53 + 05 


1 n Qyi _i_ n m 
H-U.04 it U.Ul 


AA 01 4. 1 

44. Zl it 1.00 


n f^Q 
U.Oo 




Om fl A nr K 

zuiu i\pr. 


oozy i.oi 


21.13 ± 0.04 


24.69 ± 0.05 


n qo -U n m 
+U.0Z it U.Ul 


/IK OK _|_ 1 KQ 

40. zo it i.oy 


n £.0 

U.OZ 




om n A k 

ZUIU -rt.pr. 


oozyz.oi 


21.30 ± 0.05 


23.90 ± 0.06 


1 n oa _i_ n m 

-|-U.Z4 it U.Ul 


A Q O/l 4. 1 71 
4o.y4 it 1. / 1 


n ^1 
u.01 




ZU11 rCD. lo 


OOOUO.Ol 


4.45 ± 0.04 


4.17 ± 0.05 


n 1 q _|_ n 0*3 
— U.lO It U.UO 


*i 01 -1-0 11 
O.ZO It U.ll 


n ^0 

U.OZ 




oni 1 PoK 1 
zuii ruu. 10 


oouu ( .yo 


4.07 ± 0.04 


3.87 ± 0.05 


n 1 n -L n oq 

U.1U ^ U.UO 


9 KQ _L n HQ 

z.oy u.uy 


U.O ( 




2011 Feb. 18 


55610.08 


3.85 ± 0.10 


3.98 ± 0.10 


+0.07 ± 0.07 


2.03 ± 0.07 


0.59 




2011 Feb. 20 


55612.97 


3.31 ± 0.05 


3.84 ± 0.05 


+0.30 ± 0.04 


1.41 ± 0.05 


0.58 




2011 Feb. 24 


55616.69 


2.54 ± 0.04 


2.95 ± 0.05 


+0.30 ± 0.05 


0.966 ± 0.035 


0.59 




2011 Feb. 27 


55619.73 


1.75 ± 0.10 


2.42 ± 0.08 


+0.66 ± 0.13 


0.760 ± 0.028 


0.57 




2011 Mar. 3 


55623.66 


1.32 ± 0.04 


1.64 ± 0.05 


+0.44 ± 0.09 


0.610 ± 0.023 


0.55 




2011 Mar. 7 


55627.10 


1.11 ± 0.05 


1.26 ± 0.10 


+0.25 ± 0.19 


0.500 ± 0.019 


0.55 




2011 Mar. 9 


55629.66 


1.13 ± 0.05 


1.38 ± 0.08 


+0.41 ± 0.15 


0.460 ± 0.018 


0.55 




2011 Mar. 20 


55640.04 


0.63 ± 0.03 


0.74 ± 0.04 


+0.33 ± 0.15 


0.240 ± 0.042 


0.48 




2011 Mar. 22 


55641.95 


0.62 ± 0.05 


0.76 ± 0.05 


+0.41 ± 0.19 


0.16 ± 0.04 


0.55 




2011 Apr. 15 


55666.89 


0.30 ± 0.02 


0.40 ± 0.03 


+0.58 ± 0.20 


(3.1 ± 1.0)xl0~ 2 


N.A. 


Swift/XRT 


2011 Apr. 27 


55678.84 


0.88 ± 0.10 


0.39 ± 0.10 


-1.03 ± 0.54 


(1.4 ± 0.4)xl0" 2 


N.A. 


Swift/XRT 


2011 May 12 


55693.89 


Combined upper 


limits: < 0.060 


N.A. 


(1.78 ± 0.10) xlO" 3 


N.A. 


Chandra/ ACIS 


2011 Sep. 9 


55813.33 


Combined upper 


limits: < 0.054 


N.A. 


(2.2 ± 1.0)xl0" 3 


N.A. 


Swift/XRT 


2012 Jun. 24 


56102.75 


0.30 ± 0.04 


0.35 ± 0.05 


+0.31 ± 0.39 


(5.4 ± 1.8)xl0- 2 


N.A. 


Swift/XRT 


NOTE: C03: see 


alsolCorbel 


et al.N2003h. N05: 


see alsolNowak et alj ll2005f). The calendar dates correspond to the the mid-point 



of the radio observations. Upper limits are given at a 3<r confidence level. The Galactic ridge emission has been subtracted from 
the quoted unabsorbed 3—9 keV X-ray flux. See end of section ^. 2l for the definition of the X-ray hardness. GBM: X-ray flux from 
RXTE/PCA Galactic Bulge Monitoring of Craig Markwardt. Swift/XRT: X-ray flux from Swift/XRT. Chandra/ACIS: X-ray 
flux from Chandra/ACIS. S.x.x: See complementary information in Section x.x. Na: Data points from original C03 sample, but 
these two points with optically thin radio spectra are representative of the re-ignition of the compact jets during the soft to hard 
state transition (e.g. see section [3.41 1. 



ilar X-ray flux (an unabsorbed 3-9 keV flux of ~ 2 x 10 
erg s _1 cm -2 with an error of the order of ~ 10%). 



No detection was achieved in the combined 9 GHz 
dataset with an rms noise level of 17 /xjy, whereas the com- 
bined 5.5 GHz provided a detection of GX 339—4 at the 
level of 73 ± 16 /iJy (noise level in the image of 10 /iJy) . All 
together, the 5.5 and 9 GHz data (but with a significantly 
higher noise level at 9 GHz) gives a detection of 45 ± 16 
/xjy. The choice of the final radio flux in this very faint state 
depends of the assumed spectral shape: 45 ± 16 (j,Jy if the 
spectrum is flat or 73 ± 16 pJy at 5.5 GHz in case of steep 
spectrum (which will translate to 54 /ijy at 9 GHz for a -0.6 
spectral index). For the rest of the paper, we will assume a 



Dunn et al. l2008t);lDroulans et al. 


hOld); Wuet al. d2010Tl; 


Motta et all (l201lT); IStiele et al.1 ( 


2011): ICadolle Bel et al.1 


(201lf); Rahoui et al.1 120121): Buxton et al. (20121). 



3.2 GX 339 4 close to quiescence 

Following the CABB upgrade (April 2009) of the ATCA, 
three of our radio observations have been conducted when 
GX 339—4 was in a very faint state of X-ray emission (2009 
October 30, 2011 May 12 and 2011 September 9), prior and 
after the 2010/11 outburst. Although these radio observa- 
tions did not result in a positive detection when consid- 
ered individually (Table 1), we decided to combine them 
together as the simultaneous Swift observations gave a sim- 
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Figure 4. 9 GHz radio emission from GX 339—4 in the hard 
state versus the un-absorbed 3-9 keV flux for the new data- 
points (empty green hexagons) and the ol d data-points (fille d 
red hexagons) from the original correlation ijCorbel et al,|[2003t ). 
The dashed line illustrates the original radi o/X-ray correlation 
in GX 339-4 fro m the 1997-1999 period ijCorbel et al.1 120031: 
iNowak et al]|2005l) . The full line corresponds to a fit to the new 
whole sample with a function of the form F Rad oc (with 
b = 0.62 ± 0.01). The blue point corresponds to the level of 
GX 339—4 close to quiescence (see Section l3~2ll . 



level of 45 ± 16 /iJy, which is consistent with the expected 
spectral shape of compact jets. 

Whether or not it represents the true level of 
GX 339—4 in quiescence, it is however interesting to 
note that its X-ray emission has been hovering around 
the same value for all X-ray observations c onducted be- 
tween outbursts during the la s t deca d e (e.g. Corbel et alj 
20031; iGaUo. Fender fc Corbell 120031 ; ICorbel et al.l 120051 ; 
Yen fc Konell2009l). The faintest upper limi t recently re- 



ported bv lMaccarone. Russell fc Lewis! (|2012l ) is only a fac- 
tor ~ 2.5 lower than the above X-ray flux, meaning that 
it should at least represent the level of X-ray emission, to 
within a factor 2, between the major outbursts. It may also 
be close to the level in quiescence, indicating that GX 339—4 
may be at the upper end of the quiescent BH luminosities 
(|Corbel. Tomsick fc Kaaretll2006h . 



3.3 The overall flux correlations 

3.3.1 The Radio /X-ray flux correlation 

The new radio and X-ray measurements of GX 339—4, taken 
during the la st ten years, allow u s to complement the origi- 
nal study of ICorbel et al.l (|2003l ) by adding measurements 
from seven additional outbursts, including four new out- 
bursts that encompass a full X-ray states cycle (i.e. a full 
path in a given HID). In addition, it also samples the initial 
and final hard states for a given outburst. These measure- 
ments are the most complete data-set available to study the 
evolution of the radio /X-ray correlation o f a BH over sev- 
eral outbursts (see also lCoriat et al.ll2011bl for H 1743-322). 
This new dataset represen ts an increase of the original sam- 
ple of ICorbel et al.1 (120031 ) by almost a factor 7 (88 obser- 
vations instead of 13), plus 12 observations close to state 
transitions. 



Figure 5. Same as Fig. |4]but with the data-points separated by 
the type of hard state within outburst: rising hard state (filled 
red hexagons), decaying hard state (empty green hexagons) and 
outburst in persistent hard state only (filled blue hexagons). The 
full line corresponds to the fit to the new whole sample. 



The whole sample is represented in Fig. [4] where we 
separate the new data- points (in gr e en) fr om the original 
sample (red points) of ICorbel et al.l (|2003l ). In Fig. [5] we 
divide the full sample according to the rise and the decay of 
the outbursts, as well as highlighting those outbursts when 
GX 339—4 was found in the hard state only (see also Fig. 
[H]for the data associated to each individual outburst). 

The new data-points lie very close to the original cor- 
relation, illustrated by the red points in Fig. [4] and the up- 
dated fit (dashed line) to the old data (taking into account 
the upgrade of the RXTE/PCA softwares and calibration). 
In addition, they extend the correlation further to lower and 
higher fluxes. The new sample covers almost five decades in 
X-ray flux and confirms that the correlation between radio 
and X-ray fluxes in the hard state of GX 339—4 is indeed 
very strong. The correlation is remarkably stable over a pe- 
riod of 15 years, even though several full outbursts (sepa- 
rated by off states) occurred during this period, which might 
have affected the structure/geometry of the accretion flow. 

The larger sample allow us to characterize more pre- 
cisely the functional form of this correlation. If we note, 
J*rad, the radio flux density at 9 (or 8.6) GHz (in mjy) and, 
Fx, the unabsorbed 3-9 keV X-ray flux (in unit of 10 -10 
erg s _1 cm -2 ); a fit (not using the upper limits) to all hard 



state data-points with a function of the form: 



Frcui — a, Fx 



(1) 



led to a = 1.85 ± 0.02 and 6 = 0.62 ± 0.01 (to compare with 
a = 1.99 ± 0.04 and b = 0.66 ± 0.02 for the 1997-99 period 
taking into account the RXTE using up-to-date instrument 
response files). All reported errors for the fitted functions 
are provided at the 90 % confidence level. The detection of 
GX 339—4 at low flux (Section I3.2f) . represented by a blue 
point in Fig. 3J is also consistent with this function, mean- 
ing that the global correlation could be maintained down to 
quiescence if the level we measured is really representative of 
this low flux state. The derived index of the correlation for 
GX 339-4 is fully co nsistent with the value (0 . 63 ± 0.03) re- 
obtained recently by IGaUo. Miller fc Fender! (|2012l ) for the 
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Figure 6. 9 GHz radio emission (green points) and H band (1.6 
/im) near-infrared emission (red points from Coriat et al. 2009) 
from GX 339—4 versus the un-absorbed 3-9 keV flux. The solid 
lines represent the fitting function according to the sample. For 
the infrared data, the dashed line illustrates the break at low 
X-ray flux. 



The total magnitudes have been converted to flux 
units using a total extinction of Ay = 3.7 mag. (see 
ICoriat et all 120091 ; iBuxton et all |2012| , for more details). 
With a radio to infrared spectral index that is usually pos- 
itive, it is not surprising to find brighter infrared emission 
than radio emission for a given X-ray flux (we note that 
the accretion disc may eventually contribute slightly to the 
near-infrared emis sion at low flux). As demonstrated by 
ICoriat et ail (|2009h . a broken power-law was necessary to fit 
the near-infrared/X-ray sample (x 2 = 584 for 130 degrees of 
freedom), and here the radio/X-ray sample was fitted with a 
single power-law (x 2 = 6001 for 94 degrees of freedom). We 
can not constrain the presence of a break in the radio/X-ray 
flux correlation due to the low number of radio detections 
at low flux density. The resulting reduced x 2 further indi- 
cate that intrinsic variability is also present along the fitted 
functions. The variance of the residuals (data over model) 
is 0.026 for the near-infrared/X-ray sample and 0.14 for the 
radio/X-ray sample. This highlights (see also Fig. [6} a larger 
dispersion for the radio data along the fitted function com- 
pared to the near-infrared sample. 



upper track of the correlation (their sample does not include 
our new GX 339—4 data). 

By looking at Fig. [4] which compares the original and 
new measurements, we note a significant dispersion (beside 
measurement errors) of the data-points around the fitted 
function. This dispersion is also well illustrated by Fig. [5] 
where data-points are separated according to the phase (rise 
or decay) of an outburst We recall that data-points from 
the original correlation mostly originated from the decay of 
the 1998-99 outburst. The new data allow us to sample the 
rising hard state (not observed before with ATCA), that 
is known t o reach a higher luminosi ty than the decaying 
hard state (|Maccarone fc Coppill20o"3l ). As expected, higher 
levels of radio and X-ray emission are observed during 
the rising hard state. The higher dispersion above ~10 -10 
erg s _1 cm -2 is partly related to differences between the 
rise and decay of outbursts. Indeed, the decaying hard state 
tends to be more radio bright for a given X-ray luminosity 
(or more X-ray faint for a given radio luminosity). This will 
be further discussed in Section 13.41 on several individual 
outbursts. 



3.3.2 Comparison with the OIR/X-ray flux correlation 

The radio emission from GX 339—4 in the ha rd state orig- 
inate s from the self-absorbed compact jets (|Corbel et alj 
l2000h . With an inverted radio spectrum (i.e. a positive spec- 
tral index, a, and a definition of the flux density, S u , as 
S v oc v a ), the synchrotron emis sion from the compact 
jets also reaches the infrared range (ICorbel fc Fenderl 20021 ; 
iHoman et ai"1l2005l ; iRussell et al1l2006l ; iGandhi et al.ll201ll ). 
implying that a strong correlation between the near-infrared 
and X-ray fluxes should also exist. This has indeed been 
demonstrated by ICoriat et al.1 d2009l) for the 2002-2007 pe- 
riod (see also lHoman et al.ll2005l ). 

We plot in Fig. [6] our new radio/X-ray sa mple together 
with the entire infrared/X-ray sample from ICoriat et al.l 



3.4 Focus on some individual outbursts 

We now turn to a more detailed discussion of three outbursts 
that have been well covered in radio and in particular the 
early phase of the decaying hard state when the compact 
jets are building up. It is important to understand the na- 
ture of the larger dispersion (not of statistical origin) of the 
radio data, especially in light of possible difference in the 
nor malisation between the rising and decaying hard states 
(e.g. IRussell et al.ll2007l , for XTE J1550-564). 



3.4.1 The 2004-2005 outburst 

In early 2004, G X 339—4 began a new outburst 
(|Buxton et al.ll2004l ), which ended around mid-2005. The ra- 
dio coverage of the hard states was much better than during 
the previous outbursts with 8 observations during the rise 
and 9 observations during the decay phase. The initial hard 
state was covered from 200 4 February 3, i.e. b efore the X-ray 
detection of the outburst (|Smith et a.l.ll2004h and ended on 
March 19 (close to the top of the right branch of the HID). It 
then spent almost five months in this part of the HID before 
starti ng a transition to a soft e r state around mid-August 

2004 (|Buxton fc Bailvnl l2004al; iHomanl |2004| ; iBelloni et all 
l2006h . According to lBuxton fc Bailvnl (|2005h . GX 339-4 un- 
derwent a transition back to the hard state around mid- April 
2005. When the first ATCA observation was conducted on 

2005 April 21, the radio emission from GX 339—4 was al- 
ready decaying with an inverted spectrum (spectral index of 
0.14 ± 0.03), suggesting the full reactivation of the compact 
jets. 

Figure [7] (a) illustrates the location of the radio and X- 
ray measurements for the 2004/05 outburst compared to the 
fit for the whole 1997-2012 sample. In early February 2004, 
at the faintest level of radio/X-ray emission, GX 339—4 was 
close to the global correlation. As the outburst evolved dur- 
ing the rising hard state, GX 339—4 moved significantly be- 
low this track. For the decay, a slightly different behaviour 
can be noticed once GX 339—4 moved back to the hard 
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state. In the initial part of the 2005 decay, at the brightest 
level of radio/X-ray emission, GX 339—4 was slightly above 
the global correlation. The difference between the rise and 
decay tends to lower at fainter X-ray fluxes. For this de- 
cay, the radio spectra are always consistent with optically 
thick synchrotron emission, which characterises the compact 
jets. The correlations are therefore not affected by potential 
"relic" optically thin synchrotron emission from previously 
ejected plasma that c ould interact with the ambient medium 
(Corbel et al.ll2010bl ). 

If we compare the brightest measurements during the 
2004 rise and the 2005 decay of the outburst, we find that 
the decaying hard state is significantly more radio bright for 
a given X-ray flux (or less X-ray bright for a given radio 
flux density) compared to the rising hard state. We find a 
difference of a factor ~ 2.25 in normalisation if we fix the 
slope to the value obtained for the full sample (see dashed 
lines in Fig.0. The difference is particularly important if we 
consider the soft to hard state transition, once the compact 
jets are being rebuilt. This slow evolution of the correlation 
is clearly one of the reasons for the scatter in the radio/X- 
ray correlation. This behaviour will be discussed in more 
detail in Section [4] 

34.2 The 2007 outburst 

Compared to the previous outbursts, the one that oc- 
curred in 2007 was rat her short. Indeed, with a reactiva- 
tion in November 2006 (|Swank et al.ll2006l and Section [23} , 
the transition to the soft stat e occurred in March 2007 
' Motta, Belloni fc Homanll2009f). with a transition back to 



• Rise 2004 
O Decay 2005 



the hard state after mid-May ( Kalemci et alj|2007t ). The X 
ray flux decreased up to mid- June, when a new hard state 
reflare was observed (see Fig. [3}. It was ass ociated with a 
re-br i ghtening in radio, infrared and X-ray jTomsick et alj 
l200Sl ; lCoriat et al.ll2009h . 

One radio observation was conducted at very faint 
fluxes close to the detectable onset of the outburst (Sec- 
tion [2131), whereas two other radio observations occurred in 
the brightest portion of the rising hard state (on the same 
day - February 4 - but separated by 12 hours), i.e. close 
to the top of the right branch in the HID. These three ob- 
servations are plotted (filled red hexagon) in Fig. [3(b). We 
can simply note that they are consistent with the global 
radio/X-ray correlation. One may possibly see, at high X- 
ray flux, a potentially rapid increase in the radio emission 
from GX 339-4 (see also the 2002 and 2010 rises in Fig. 
that could possibly be related to the onset of a radio flare 
(e.g. iGallo et alJl2004h . Howe ver, those flares usu ally occur 
at much lower X-ray hardness (I Fender et al I I 19991 ) (it is also 
likely the case in GX 339—4, e.g. see section I3.4.3[l and is 
probably not a valid explanation for the observed increase 
here. 

The decay in 2007 has also been well sampled in term of 
radio observations with the ATCA. Indeed, for the first time, 
the onset of the compacts jets has been observed during the 
transition from the intermediate to the hard state (see also 
ICorbel et al.|[2000l and the 2011 decay). The radio spectra 
are initially consistent with optically thin synchrotron. But 
as the radio emission rapidly increases (see black hexagones 
in Fig. [7] b) , it becomes and stays optically thick. As the 
focus of this paper is the long term study of the radio X-ray 
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Figure 7. 9 (or 8.6) GHz radio flux density from GX 339—4 ver- 
sus the un-absorbed 3-9 keV flux for the X-ray rise (filled symbols) 
and the X-ray decay (empty symbols) for the three outbursts of 
GX 339—4 with good radio coverage during both the X-ray rising 
and decaying hard states. The epoch of outburst are highlighted 
in the top left corner of each panel. The time of the reflare in 
2007 is also indicated in the middle panel. The solid line illus- 
trates the fit to the whole 1997-2012 sample with a function of 
the form Fpt a[ j <x (as discussed in Section I3.3.1H . The dashed 
lines illustrate the variation in normalisation of the fitted function 
needed to accommodate the observed flux variations between the 
rise and the decay with the slope fixed to the value obtained for 
the whole sample. The black hexagons and squares in the mid- 
dle and bottom panels illustrate the turn-off (filled) and turn-on 
(empty) of the compact jets. Note that the axis ranges are not 
identical in all figures. 
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correlation in GX 339—4, we will discuss this aspect together 
with the behaviour at different wavelengths in a forthcoming 
paper. We note that, during the formation of the compact 
jets, the 3-9 keV X-ray emission from GX 339—4 steadily 
declined as usual in the hard state. Once the radio emission 
is fully consistent with the compact jets, the radio emission 
also gradually decreases. 

In terms of the radio/X-ray correlation, the data points 
during the decay (Fig. [7| (b)) are indicated by the open red 
hexagons, whereas the open red triangles highlight the re- 
flare during this decay. Despite only five radio observations 
for this standard decay (as the re-flare occurred soon after 
the hard state transition), we again note that for a given 
X-ray flux, GX 339—4 appears more radio bright during the 
decay compared to the rise. A change in the normalisation of 
a factor at least 1.94 would be needed if we want to move the 
lower track (associated to the rising hard state) to accommo- 
date the brightest radio observation in the decay (keeping 
the slope fixed to the value obtained for the full sample). 

Regarding the re-flare portion of the 2007 activity, we 
can simply note that all data-points are strongly correlated 
with a slope consistent with the global correlation. A ra- 
dio observation conducted on November 27 after the end of 
the 2007 re-flare (black diamond) is also indicated. It corre- 
sponds to a period of very weak residual hard state activ- 
ity from GX 339—4. It is again globally consistent with the 
overall correlation. 



3.4.3 The 2010-2011 outburst 

A new full outburst began i n 2010 as first note d by 
the MAXI X-ray monitorin g (lYamaoka et all l2010f). Af- 
ter an initial hard state l)Tomsickl l2010l : ICorbel et all 
2010ah. a transition to the soft sta te was reported by 



Belloni. Motta fc Munoz-Dariasl (|201Cl ) . The transition back 
to the hard state occur r ed around mid- Janua ry 2011 
ijMunoz-Darias et all 120111 iRussell fc Lewis! l201ll l Com- 
pared to previous outbursts, the rising hard state has been 
well covered with 20 radio observations conducted almost 
up to the transition to the soft state. The very early part of 
the rise has, however, not been sampled this time. A similar 
number of observations was conducted for the decaying hard 
state, meaning that it is the best sampled outburst with the 
most sensitive (due to CABB) ATCA radio observations. 

Inspection of the radio and X-ray measurements in 
Fig. [7] (c) shows that they are rather well correlated and 
consistent with the global correlation discussed in the pre- 
vious subsections. However, as the outburst evolves in the 
rising hard state, we note that GX 339—4 seems to get more 
radio bright (or less X-ray bright) for a given X-ray flux, 
possibly implying a steeper correlation at higher flux. This 
steepening of the correlation at higher flux (very close to 
the peak flux in hard state) may possibly be a recurrent 
phenomena, as similar behaviour seems to be present also 
in 2002 (Fig.EJ and 2007 (see Section l3A2j) . At the end of 
the hard state, the r adio emission starts to be quenched a s 
observed previously (jFender et al.ll 19991 : ICorbel et alll2000h . 
It also indicates that the major radio flare in GX 339—4 oc- 
curs at a hardness significantly below the one defined by the 
hard state. 

The transition back to the hard st ate occurred at lower 
X-ray flux defining an hysteresis cycle (|Maccarone fc Coppil 



l2003h - Similarly to 2007, the radio emission is initially opti- 
cally thin (open black squares in Fig. c). The radio emis- 
sion then becomes and remains optically thick (a signature 
of a full re-building of the compact jets) once GX 339—4 
gets close to the global radio/X-ray correlation. Again, sim- 
ilarly to the decay in 2005 and 2007, it is initially above 
the standard correlation and as the X-ray emission decays, 
GX 339—4 is observed to move back to the global correla- 
tion. In the very bright parts of the hard state (and close 
to the state transition), significant thermal emission from 
the accretion disc could contribute to the X-ray spectrum. 
If we only keep the contribution from the power-law, the 
points close to the transition (black squares in Fig. [7] c) 
move obviously towards lower X-ray fluxes, but not back to 
the global correlatio n (as it may be the case in Cyg X— 1; 
IZdziarski et allboill ). 

If we adjust a line to fit the two extremes in the radio 
and X-ray diagram (lower track for the rising hard state in 
2010 and upper track for the decaying hard state in 2011) 
but with a slope fixed to the global correlation, we need 
a factor 2.19 (identical to 2004-2005) to accommodate the 
observed variations. Again, the larger dispersion along the 
global correlation occurs around the time of compact jet 
formation. 



4 DISCUSSION 

We have presented a set of 100 quasi-simultaneous radio and 
X-ray observations of GX 339—4, including 12 observations 
close to state transitions. This constitutes the largest sample 
of measurements over a broad range of fluxes and also over 
the longest time period for a black hole in the hard state. 
In addition, this long term study of GX 339—4 allows us 
to present the evolution of the radio/X-ray flux correlations 
over several outbursts separated by off states for a given BH 
without any bias from distance uncertainties. This makes 
GX 339—4 the reference source for comparison with other 
accreting sources (BHs, neutrons stars, white dwarfs and 
AGN). The main characteristics from the new set of data 
can be summarised as follows: 

• Despite numerous outbursts of different nature (full or 
hard state only) during a 15-year period, the overall corre- 
lation between the radio and X-ray fluxes is strongly main- 
tained. It is observed over five decades in X-ray flux and 
three decades in radio flux. By fitting all measurements (ra- 
dio flux density at 9 GHz in mjy, Fl-ad, and unabsorbed 3-9 
keV X-ray flux in unit of 10 -10 erg s _1 cm -2 , Fx) with a 
function: F la a = a, F x , we obtained a — 1.85 ± 0.02 and 
6 = 0. 62 ±0.01. The faintest levels of X-ray and radio emis- 
sion, possibly GX 339—4 in quiescence, are also consistent 
with the global correlation. 

• We note a significant dispersion around the fit- 
ting function in the radio/X-ray diagram that appears 
more pronounced than for the IR/X-ray correlation (from 
ICoriat et all 120091 ). However, the dispersion of the radio 
measurements does not occurr randomly. The strongest de- 
viation seems to occur in periods when the compact jets are 
turning on, during the soft to hard state transition (see 2005, 
2007 and 2011 decaying hard states, see Fig. [7J). Additional 
deviations may exist at higher flux during the rising hard 
state. 
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Figure 8. 9 GHz radio emission from GX 339—4 versus the un-absorbed 3-9 keV flux for the whole 1997 - 2012 period. The different 
outbursts are highlighted, as well as some individ ual portions ( e.g. ri se or decay) if necessary. The lines correspond to the estimated (not 
a fit) radio flux densities based on the relation of lCoriat et alj (120091) using the observed near-infrared/X-ray correlation (valid for X-ray 
fluxes above 3 X 10~ n erg s _1 cm~ 2 , a region highlighted by the thicker lines), only assuming a radio to near-infrared spectral index, 
ctR at i-iR, ranging from 0.25 to 0.05 with a step of 0.05. The extrapolation to lower fluxes is indicated by the dotted lines (but see text 
for cautions). 



• During the period of compact jet formation in the de- 
caying hard state ( < 3% Lsdd for a 6 Mq BH at 8 kpc), we 
note that GX 339—4 seems to be more radio bright (or less 
X-ray loud) than during the corresponding rising hard state, 
resulting in a difference of normalisation of the order of 2, 
that seems constant over the outburst history of GX 339—4. 



4.1 A coupled radio, near-infrared and X-ray 
correlation 

In Corbel et al. (in prep.), we study in detail the formation 
of the compact jets during the soft to hard state transi- 
tions. We show that the compact jets first turn on at radio 
frequencies and then in the QIR (a conclusion th at has also 
recently been reached bv lMiller- Jones et al.| [2012). However, 
the timescale to reach the peak is much longer in the near- 
infrared compared to the ra dio domain (see for example also 
Fig. 5 in lCoriat et al.ll2009h . The consequence is that the ra- 
dio emission quickly decays, whereas the near-infrared emis- 
sion is still rising, when the 3-9 keV X-ray flux is steadily 
decreasing. This implies a significant evolution of the radio 
to near-infrared spectra (equivalent to major changes in the 



radio to near-infrared flux ratio, e.g. see ICoriat et aI1l2009h 
during the soft to hard state transition. 

The dispersion in the near-infrared / X-ray flux correla- 
tion is smaller than in a comparison between radio and X-ray 
emissions (Section I3.3.2|l . This is consistent with the stan- 
dard picture where the NIR emission zone is located close to 
the base of the jet and is therefore more tightly connected 
to the X-ray emitting flow than the termination of the com- 
pact jets where the radio emission arises. In that case, the 
exact location of GX 339—4 in the radio / X-ray flux dia- 
gram would then be influenced by the radio to near-infrared 
spectral index (or similarly by the radio to near-infrared flux 
ratio), i.e. the internal energetic balance of the compact jets. 

Indeed, in a generic mod el of self absorbed compact 
jets (|Blandford fc Konid[l979l ), the jet emission is charac- 
terized by a flat to inverted spectrum from the radio up 
to a turnover frequency, that is estimated to be in the mid 
to near-infrared range f o r Galactic BH s dCorbel fc Fenderl 



20021 ; ICoriat et ail 120091 ; iGandhi et ail l201ll ; iRahoui et al 



20111 ). In this optically thick regime, the radio, -Fkad, to 
near-infrared, -Fkear-m, flux densities can be written as : 
Jriad cx F Nc ar-iR !/* Ri » d - IR , with «R a d-iR the radio to near- 
infrared spectral index and v the frequency. Using the re- 
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lation of lCoriat etHI <|2009h between the near-infrared flux 
densities and the 3-9 keV X-ray fluxes, F x , of Jkear-iR = k 
Fx 2 (with 62 = 0.48 ± 0.01 for the portion of the correlation 
at high flux), we obtain the following simplified relation: 



l/Rad 
I'Near-IR 



(2) 



with fR a d and ^Noar-m the frequencies corresponding to the 
radio and near-inf rared data (k is adapted from Table 1 of 
ICoriat et aj]l2009l ). 

However, a constant radio to near-infrared spectrum 
(i.e. a_R ac j_/_R=constant) should normally lead to a depen- 
dence between radio and X-ray emission of the same kind 
as observed between infrared a nd X-ray emission, i.e. FR a( j 
tx F x 2 (with b 2 = 0.48 ± 0.01. ICoriat et al.ll2009h . This is 
not what is exactly observed as we found -pR a( j oc F x , with 
b = 0.62 ± 0.01 (Section [3XTJ. But, equation © also high- 
lights an extra dependance of the normalisation with the 
radio to near-infrared spectral index ctRad-iR, that could 
influence the exact location of GX 339—4 in the radio/X- 
ray flux diagram (beside the X-ray flux) and could affect 
the fitted correlation index b. 

As a complementary comparison, we plot in Fig. [8] the 
simultaneous radio/X-ray measurement for the various out- 
bursts of GX 339—4. The different lines (not a fit) corre- 
spond to the predicted radio fluxes based on t he measured 
infrar ed /X-ray correlation (equation [5J) from ICoriat et alj 
|2009l ) (we use the functional form with no break that is 
valid for 3-9 keV flux above ~ 3 x 10 -11 erg s _1 cm' 2 ), 
assuming a radio to near-infrared spectral index QR a d-m 
in the range of 0.05 to 0.25, which are typical values for 
the hard state. These small changes in the radio to near- 
infrared spectral index are sufficient to explain the overall 
dispersion in the radio and X-ray flux correlation, especially 
the larger dispersion above ~10 -10 erg s _1 cm~ 2 (because it 
corresponds to larger variation in QR a d-iR, see Corbel et al. 
in prep.). Such an evolution of the radio to near- infrared 
spectral index could simply be due to deviation from the 
simple Blandford-Konigl model used for compact jets (such 
a different jet g eometry and/or different electron energy dis- 
tributions, e.g. lHeinz]|2006l ; [Ka iser 2006]) or the efficiency of 
particles acceleration in the jets (Corbel et al. in prep.). 

However, one needs to be cautious with equation [2l 
which assumes that the near-infrared observing band lies 
on the optically thick part of the jet spectrum. While this 
should be fine at high flux (although not really tested), this 
may not be the case at lower flux (see ICoriat et al. 
for more details). Furthermore, the near- infrared emission 
of GX 339—4 at high flux in the hard state is dominated by 
the compact jets (e.g. the large drop in OIR emission dur - 
ing the hard to sof t state transitions. | Homan et al.l (|2005h ; 
ICoriat et all (|2009l ); iBuxton et al.l (|2012l )). but the thermal 
emission from the accretion disc may contribute at low flux 
in the total near-infrared emission from GX 339—4. There- 
fore, the thick lines in Fig. [8] should be a guide to highlight 
a possible link between radio, near-infrared and X-ray emis- 
sion, but the full understanding of the radio to near-infrared 
SED is beyond the scope of this paper. 



4.2 The "Universal" radio/X-ray flux correlation 

4.2.1 The "standard" track 

The new sample of GX 339—4 provides a new opportunity 
to compare the radio and X-ray emission of Galactic accret- 
ing binary BHs and neutron stars. In Figure |3 we report 
our new data of GX 339—4 along with the measurements 
from published studies representing a total of 24 BHCs. 
This includes the very few sources with very good radio and 
X-ray cove rage (i.e. V404 Cyg from Gallo. Fender fc Poolevl 
(|2003h andlCorbel. Koerding fc Kaaretl (|2008l ); H 1743-322 



from ICoriat et al 



(|2011bh : GRS 1915 +105 from 
iRushtonet al.l (|2010h . Swift J1753. 5-0127 fromlSoleri et all 
j2010t ). MAX! J1659 -152 from iRatti et al.l (|2012T ) and 
XTE J1752-223 from Ijonker et all (|2012T ) and Brocksopp 
et al. to be submitted) and the very few measurements 
from additional sources rep o rted in more global studie s 



(e.g. iGallo. Fender fc Poolevl ((20031); ICalvelo et al.l (120101) 



IFender, Gallo fc Russelll (|201oh ; ICoriat et al.1 (12011a!)) or 



indiv idual reports (e.g. IGR J17 177— 3656 from Paizis et al.l 
|2011| ; IGR J17091-3624 from iRodriguez et al.l [201ll V To 

avoid confusion in the figure, we do not plot the large 
group of measurements related to Cyg X— 1. Two neutron 
star systems are a ls o rep o rted for compariso n : Aql X— 1 



star systems are a ls o rep o rted tor com parison: Aql A— 1 
from iTudose et al.1 J2009h: llVIiller- Jones et al.l (|2010h and 



4U 1728-34 from lMigliari fc Fender! (120061 ). In addition, we 



also report the quie s cence levels of a sample of BH binaries 
jCalvelo et al.ll2010l : iMfller- Jones et al.ll201ll ). All X-ray lu- 
minosities have been converted to t he 1-10 keV band and we 
used the updated distances from IFender. Gallo fc Russelll 
|2010l . and references therein). For some sources, especially 
towards the Galactic bulge, the distance may be unknown 

and is assumed to be at 8 kpc. 

With the distances of G X 339-4 (IZdziarski et al.ll2004r i 
and V404 Cyg (updated in iMiller- Jones et aLlbOol h these 
two sources now share the same track in the radio/X-ray 
diagram. The updated fit to the whole GX 339—4 sample is 
represented on Figure [9] by the sol id line. The quiescent level 
of A 0620-00 (|Gallo et al.ll2006l ) is remarkably consistent 
with an extrapolation of the GX 339—4 fit to lower fluxes. 
However, the exact track of A 0620—00 during its outburst 
is unknown as no radio observatio n was conducted durin g its 
hard state at that epoch (see also lYuam Yu fc Hdl2009l ). All 
other upper limits from quiescent BHs does not contradict 
the extrapolation of the correlation defined by GX 339—4. 
The extension of the radio/X-ray correlation with index ~ 
0.6, as determined by the new data of GX 339—4, down 
to the faintest levels of emission from Galactic black hole 
binaries, implies that the standard correlation could be sta- 
ble (despite minor deviations within an outburst) over more 
than 8 decades in X-rays and 6 decades in radio luminosi- 
ties. Furthermore, no break in the correlation is observed 
at lower luminosity with the a vailable data (as e xpected in 
some theoretical model , e.g. lYuan fc Cuil 120051), but this 
could occur at even lower luminosity ( Yuan. Yu fc Hof2009l ) . 
This means that more sensitive observations (e.g. with the 
upgraded JVLA/ATCA and/or the new forthcoming SKA 
pathfinders) of quiescent BHs (or BHs in their decaying hard 
state) would be needed to better constrain the slope of the 
correlation at low luminosities, as it currently relies only on 
the quiescence level of A 0620—00, V404 Cyg and possibly 
now also GX 339-4. 
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Figure 9. Radio and X-ray (1-10 keV) luminosities for Galactic accreting binary BHs in the hard and quiescence states (see text for 
references). It illustrates the standard correlation (defined by sources such as GX 339—4 or V404 Cyg with index ~ 0.6) and the new 
correlation for the outliers (defined by e.g. H 1743—322 or Swift J1753.5— 0127 with index ~ 1.4 ). The solid line illustrates the fit to the 
whole 1997-2012 sample of GX 339—4 (as discussed in Section |3. 3. H with an extrapolation to th e quiescence state of BHs. The dashed 
line corresponds to the fit to the data for H 1743—322, one of the representatives for the outliers ijCoriat et alj|201 lbl l ■ Upper limits are 
plotted at the 3cr confidence level. For a few sources, their distances are unknown. 



4-2.2 The "outliers" track 

A growing number of sources (originally XTE J1650— 500 
(ICorbel et all 12004) . IGR J1749 7— 2821 ([Rodriguez et all 



x 10 36 erg s" 1 the outliers seem to leave the H 1743-322 
track and join the standard correlation around 10 35 erg s _1 . 



2007ft, Sw i ft J 1753. 5-0127 (|Cadolle Bel et all l2007t 



Soleri et all l201Ch ) are now detected significantly below 
the standard correlation. The number of these outliers is 
rapidly increasing (see Figure |5). At a 1-10 keV luminosity 
above 4 x 10 36 erg s _1 , the group of outliers seems to be 
consistent with the track followe d by H 1743—322 w ith a 
steeper correlation index of ~ 1.4 (|Coriat et al.ll2011bft . and 
also quite similar to the few reported neutron star systems 
l|Migliari fc Fenderl 12006ft . In any-case, this larger dataset 
confirms the existence of two tracks in the radio/X-ray flux 
correlation at higher X-ray fl ux, but still in the hard stat e 
(recently statically proved bv lGalloT Miller fc Fendeill2012l) . 
It is unclear for the group of outliers if the normalisation 
of the correlation is similar or not, as some dispersion is 
also obser ved along the H 1743 —322 track. Furthermore, as 
found by ICoriat et all (|2011bl lah and reinforced here with 
the increased numbe r of sources (e.g. MAXI J1659— 152 
and XTE J1752-223 iRatti et al.ll2012l ; Ijonker et ail 120121 
Brocksopp et al. to be submitted), at luminosity below 4 



The spin of the BH (|Fender. Gallo fc Russelll2010l) . as 
well as the binar y parameters (orbital period, accretion disc 
size, inclination, ISoleri fc Fendej|2MTl ) do not seem to play 
a role in defining the level of jet power. By using a jet toy- 
model whli_aii_jncreasin^_bulk Lorentz factor above 10 -3 
ISoleri fc Fender] (|201lD were able to reproduce the 
larger scatter for the brighter hard state (the zone of the 
outliers). It is, however, clear from Figure [9]that the outliers 
are a separate population from the sources lying on the stan- 
dard correlation (see also lGallo. Miller fc Fenderj|2012T ). The 
radio/X-ray diagram is not uniformly populated, as would 
have been expected if the two populations were related to 
a uniform distribution of inclination angle (and henceforth 
Lorentz factor). The existence of a critical magnetic field 
(ICasella fc Pe has also been invoked to explain the 

lower radio luminosity for the outliers, but then it also needs 
to explain the existence of two populations as well as the 
steeper correlation index. In order to find an explanation 
of the outliers by a lower radio luminosity (or jet power), 
it may then be necessary to invoke a different coupling be- 
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tween jet power and mass accretion rate (see ICoriat et alj 
l2011bl for more details). 

Instead of characterising the ou tliers by fainter jet 
power (as in iSoleri fc Fended 1201 ll ). one may instead 
consider these sources as being more X-ray bright than 
the ones from the standard correlation, meaning that the 
understanding of the correlations would be related to the 
origin of the X-ray emission process(es). We would then be 
left with two groups of sources above 4 x 10 36 erg s _1 (one 
with an index of ~ 1.4 and the other one with the typical 
index of ~ 0.6). It seems apparent that only th e ~ 0.6 track 
exists a t lower radio and X-ray luminosities. ICoriat et al.l 
l|2011bl . and references therein) have extensively discussed 
the nature of these two groups. In that framework, the 
outliers would be consistent with the presence of radiatively 
efficient BHs in the hard state (hot accretion flows or 
accretion disc corona). The standard correlation would still 
be related to a radiatively inefficient accretion flow (such as 
X-ray jets or ADAF models). The difference between the 
two populations could po ssibly be related t o variation in 
the a viscosity parameter (|Xie fc Yuanll2012T ). 



4-2.3 Caveats on constraining the slope of the correlation 
using portion of one individual outburst 

T he global resu l ts prese nted in this study for GX 339—4 or 
in ICoriat et ail (|2011rJ ) for H 1743-322 allow us to char- 
acterise the evolution of the radio and X-ray fluxes correla- 
tion for a given source along very different outbursts. While 
the global correlation is maintained with a specific slope (~ 
0.6 for GX 339-4 or ~ 1.4 for H 1743-322), we note that 
temporary deviations do exist. This happens for GX 339—4 
several times, mostly in the decaying hard state once the 
compact jets have been fully rebuilt (see 2005, 2007 and 
2011 for the decay), but also possibly even at higher fluxes 
(although it is not so clear at present, see end of rises in 2002, 
2007 and 2010). During these limited portions of outburst, 
correlation indices in the to 2 range can temporarily be 
observed. It also happens to H 1743—322 whi ch shows two 
track s and a transition (w ith a flatter index, Ijonker et al.l 
l20ld : ICoriat et all l201ibl ) between the two tracks. Addi- 
tionally, we note that other sources have also been ob- 
served with different corre l ation indices (e.g. GRO J1655— 40 
IShaposhnikov et al. I l2007l ; IXue fc Cuil 120071 ) . This concerns 
Galactic BHs (and probably also accreting neutron star bi- 
naries) . 

The observed variations of the correlation index in 
GX 339—4 (see Fig. [TJl occur over a 3-9 keV flux range 
of roughly a decade. Several reasons could be invoked to 
explain temporary variations of the correlation index, such 
as changes in the radio spectral index or the radio to near 
infrared spectral index (as discussed in section 4.1), an in- 
crease of the importance of electrons cooling in the jets, the 
presence of an additional X-ray emitting component, etc. In 
order to properly constrain the global correlation index, the 
3-9 keV (or a similar X-ray band) flux should be sampled 
over more than a decade in X-ray flux, and probably ob- 
servations over two decades in X-ray flux should provide a 
reasonable es timate of this inde x (e.g. see the example of 
H 1743-322. rCoriat et alj|2011bl ). However, if the source is 
an outlier (section l4.2,2[) , special care should be taken (with 



regular monitoring) in order to separate the different tracks 
(mixing the 1.4 track with observations taken during the 
transition to the standard 0.6 track could lead to a different 
index) . 

As the c orrelation has been ext ended 
to AGNs (iMerloni Heinz fc di Matteol 120031 ; 
iFalcke. Kording fc Markofj 120041 '). it also means that 
very different correlation indices could potentially be 
mea sured in superma ssive BHs with very limited datasets 
fe.g. lKing et~aLll201lh . Therefore, one has to be very careful 
before generalising the slope of the correlation from a very 
specific portion of an outburst. Only very long studies 
sampling large radio and X-ray flux variations of a given 
source (or a large sample of different sources with known 
distance) could reveal the intrinsic empirical relations 
necessary to understand the coupling between accretion 
and ejection mechanisms in accreting BHs. 



4.3 On the nature of "hysteresis" in OIR 

Durin g the 2000 outburst of XTE J1550-564, iRussell et all 
(|2007h plotted the infrared versus X-ray emission and re- 
ported quasi-parallel tracks for the rising and decaying hard 
states. This phenomenon, that they called "hysteresis", has 
been reported to date only in XTE J1550— 564. Whereas, 
the difference in normalisation betwee n the rise and deca y 
is of a facto r 5 in XTE J1550-564 (jRussell et al.ll2007f ). 
ICoriat et al.l (|2009h already pointed out that no difference 
could be seen in infrared between the rise and decay of 
GX 339—4, despite an intense monitoring over 5 years. 
However, as discussed in Section [3.41 we found some differ- 
ences in radio between the decaying and rising hard states. 
This occurs roughly during period of formation and destruc- 
tion of the compact jet (similarly to XTE J1550— 564). Al- 
though, we clearly do not have parallel tracks in the case 
of GX 339—4, we possibly observed changes up to a fac- 
tor 2 in term of variation of the normalisation of the ra- 
dio track. Similarly to XTE J1550— 564, the jet radio emis- 
sion in GX 339—4 is sometimes brighter during the de- 
cay compared to the rise for a given X-ray luminosity. The 
radio observations conducted during the 2000 outburst of 
XTE J1550-564 (|Corbel et alj|200l| ) are too sparse to be 
compared with the infrared data. 

To explain th e "hysteresis" in XTE J 1550 -564, 
IRussell et all <|2007h considered changes in the radiative ef- 
ficiency of the accretion flow, variations in the viscosity pa- 
rameter a or modification in the jets properties. Amongst 
these different scenarios, one needs to find one that could 
explain the limited "hysteresis" in GX 339—4 (only a factor 
2 in radio and no difference in near- infrared). In GX 339—4 
and XTE J1550— 564, the radio and infrared emissions in the 
hard states originate from the compact jets. The fact that 
no difference is observed in the near-infrared for GX 339—4 
possibly implies that the origin of this "hysteresis" is not 
solely related to the X-ray emission process(es), but rather 
to the compact jets. 

As discussed in Section T4. II the radio, near-infrared and 
X-ray emissions of GX 339—4 are strongly coupled. It is 
then possible that this "hysteresis" is only determined by 
the spectral energy distribution of the compact jets and 
therefore traces the jet behaviour (rather than the accre- 
tion flow). The level of "hysteresis" and the involved spec- 
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tral domain could possibly be related to the evolution of the 
break frequency where the jets becomes optically thin, but 
also to the evolution of radiative efficiency along the jets, 
the interplay between the jets and the corona and/or the 
efficiency of particles acceleration. Alterna tively the mag- 
netization of the disk ()Petrucci et al.l[2008l ) could possibly 
modify the broadband spectral energy distribution of the 
jets. 



5 CONCLUSIONS 

We reported a series of radio and X-ray observations of the 
recurrent BH GX 339—4 during the past 15 years. This 
new set of observations samples almost all luminosity levels 
(nearly 5 decades in X-ray flux), possibly down to quies- 
cence, of a BH along several different outbursts. GX 339—4 
can therefore be considered a reference for comparison with 
other accreting sources without any bias from distance un- 
certainties. 

Despite numerous outbursts of different nature, the 
overall radio and X-ray emissions display a very strong 
non-linear correlation with a coupling of the form Lx oc 
^Rad ±0 01 - This very well constrained correlation index is 
consistent with previous studies, that highlighted an index 
in the range of 0.5 to 0.7. The same correlation seems to 
be maintained from the brightest hard states to the low- 
est detected luminosity level, that could be consistent with 
the quiescence state of GX 339—4. With the reported long- 
term near-infrared /X-ray fluxes correlation in GX 339—4 
jCoriat et al.ll2009h , we further demonstrated a coupled cor- 
relation between these three frequency ranges. The level of 
radio emission could then be tied to the near-infrared emis- 
sion by the evolution of the jet broad band SED (assuming 
that the near-infrared data lies on the optically thick portion 
of the compact jet spectrum). 

We also highlighted periods of significant and higher 
dispersion along the fitted function in the radio/X-ray sam- 
ple, when compared to the near-infrared/X-ray data. This 
appeared more pronounced in periods of formation and de- 
struction of the self absorbed compact jets. For a given X- 
ray flux, the radio emission could be higher (up to a fac- 
tor 2) during the decay compare to the rise of the out- 
burst, whereas no difference was found in near-infrared 
(jCoriat et al.ll2009h . By compari ng this behaviour w ith a 
similar one in XTE J1550-564 (|Russell et al.l 120071 '). this 
"hysteresis" behaviour seems to be related to the properties 
of the compact jets. 

We incorporated our new data of GX 339—4 in a more 
global study that includes a set of 24 BHCs. We observed 
that the "standard" track, that is now well constrained by 
GX 339—4, is closely shared by a subset of these BHCs (e.g. 
V404 Cyg). The quiescent BH A 0620-00 is the faintest 
reported source ( more than 2 orde rs of magnitude fainter 
than GX 339-4, iGallo et al1l2006r i and it obeys the stan- 
dard track that is defined with our large GX 339—4 sample. 
Furthermore, we collected a growing number of sources (the 
majority in fact) that appear significantly below the "stan- 
dard" track (see also IGallo. Miller fc Fended 12012] ). those 
"outliers" may possibly share the behavior reported for 
H 1743—322, meaning that they could join the standard cor- 
relation at low and/or high X-ray fluxes. We emphasize that 



special care should be taken to constrain the index of the 
correlation in very limited sample, as this index can only be 
securely obtained with more global studies such as the one 
we conducted with GX 339—4. The location of a source in a 
given track of the radio/X-ray diagram could be related to 
the radiative efficiency of the inner accretion flow or a dif- 
ferent coupling betw een jet power and mass accretion rate 
|Coriat et al.ll2011bl ) . While the standard track has been ex- 
tended to supermassive BHs by the so-called "fundamenta l 
plane of BHs activity" jMerloni. Heinz fe di Matted 120031 ; 



iFalcke. Kording fc Markofi 2004 ). the existence a popula- 
tion similar to the "outliers" within AGNs would be a strong 
support for scale invariance of the jet- accretion coupling in 
accre ting black holes (see recent work of lBroderick fc Fended 
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